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Abstract. We review theoretical models of the early stages of star for- 
mation, in which gravitational collapse is strongly regulated by magnetic 
fields and the associated process of ambipolar diffusion. We discuss re- 

OO '. suits of numerical simulations and analytical studies of core formation 

CN ■ and collapse, which can be directly tested against observation. We also 

focus on recent data which are relevant to this theory of star formation, 
such as: observations of extended infall in protostellar cores, estimates of 

^^ , evolutionary timescales ~ 1 Myr for cores, measured mass-to-flux ratios 

of cores, and the relative alignment of polarization vectors with appar- 
ent cloud elongation. It is shown that in all of these areas, the data 
remain compatible with magnetic collapse models which lie within the 
observationally allowable range of parameter space. Other areas of in- 

*—* | terest (protostellar accretion rates and the presence of core edges) and 

Ci ■ issues that remain unresolved or under study (the role of non-thermal 

motions and cluster formation) are also discussed. Moreover, we high- 
light some differences between our model predictions and those of highly 
turbulent star formation models, and discuss how these differences can 
j^ ■ be distinguished observationally. 

"Never multiply explanations or make them more complicated than 
necessary. An explanation should be as simple and direct as possi- 
ble." 

"No more things should be presumed than is necessary." 

— William of Occam (attributed to), English monk, philosopher 

Occam's Razor (c. 1285-1349). 

1. Introduction: The Need for a Magnetic Theory of Star Formation 

A robust physical theory is one that explains observed phenomena and obeys 
the laws of physics. And, in the sense of Occam described above, it should be as 
parsimonious as possible with the number of assumptions used in its formulation. 
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There are a number of lines of evidence that strongly suggest that mag- 
netic fields play a crucial role in a comprehensive theory of star formation. For 
instance, observations of magnetic field strengths by OH Zeeman splitting in 
dense gas (e.g., Crutcher 1999) allow a determination of the mass-to- flux ratio 
M/&b in molecular clouds. When allowances are made for the random ori- 
entations of clouds and magnetic fields (e.g., see Table 1 of Shu et al. 1999), 
the observed values for M/&b in several clouds are found to be a factor ~ 2 
below or above the critical value for gravitational collapse, (M/^b)ctH (e.g., 
Mouschovias & Spitzer 1976). That the observed values for these clouds are 
so close to the critical value indicates that magnetic fields must be instrumen- 
tal in their support against self-gravitational collapse. Additionally, Crutcher's 
data set also reveals that the magnetic field strength B scales with the gas den- 
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sity p n as B oc p n , consistent with the predictions from models of equilibrium 
magnetic clouds (Mouschovias 1976) and ambipolar diffusion models of mag- 
netically supercritical core formation and collapse in subcritical clouds (Fiedler 
&; Mouschovias 1993; Ciolek & Mouschovias 1994, 1995; Basu & Mouschovias 
1994, 1995a, b). 

Another inextricable link between magnetic fields and star formation is the 
angular momentum of stars and protostellar disks. It has long been known that 
if absolute conservation of angular momentum were to hold true throughout 
the evolution of a cloud, from its formation to protostellar collapse, then star 
formation would be halted by centrifugal forces (e.g., Mestel &: Spitzer 1956). 
The solution to this conundrum is magnetic braking of clouds and cores, which 
allows the efficient transfer of angular momentum via rotational Alfven waves 
along the magnetic field, from a cloud (e.g., Mouschovias & Paleologou 1979, 
1980) or core (e.g., Basu &; Mouschovias 1994, 1995a, b) to an external medium. 
Measurements do indeed find that rotational velocities in clouds are orders of 
magnitude smaller than that which would occur under strict conservation of 
angular momentum (e.g., Goldsmith & Arquilla 1985). Observations of near 
alignment of rotation axes with polarimetrically determined magnetic field di- 
rection in clouds (e.g., Heyer et al. 1987) and starless Bok globules (e.g., Kane & 
Clemens 1995) also confirm the expected correspondence from magnetic braking 
models. Another prediction of these models is retrograde rotation of clouds and 
fragments (Mouschovias & Paleologou 1979), which has also been seen to occur 
(Clark & Johnson 1981; Young et al. 1981). 

There are various other pieces of support for the importance of magnetic 
fields in star formation, such as the existence of the magnetic flux problem of 
star formation (e.g., Mestel & Spitzer 1956), evidence that clouds are in near 
virial equilibrium, with linewidths that are related to the Alfven speed (Myers 
&; Goodman 1988; Mouschovias & Psaltis 1995), and the potential importance 
of the Parker instability in forming molecular cloud complexes (Mouschovias, 
Shu, & Woodward 1974). 

As argued above, all of these various observational facts or threads of rea- 
soning can be explained and tied logically together through a single unifying 
theory of star formation in magnetically supported clouds. In the next section 
we summarize results of ambipolar diffusion models of core and star formation. 
More extensive treatments of the history and the relevant physics of the magnetic 
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models of star formation can be found in the detailed reviews by Mouschovias 
(1987, 1996) and Mouschovias & Ciolek (1999). 

2. Key Features of Ambipolar Diffusion and Core Formation 

Ambipolar diffusion, the drift of neutral matter with respect to plasma and 
magnetic field, was originally suggested by Mestel & Spitzer (1956) as a way 
that magnetic flux can leak out of interstellar clouds and induce collapse and 
fragmentation. A fundamental reevaluation of this concept was introduced by 
Mouschovias (1976, 1977, 1978), who pointed out that the main effect of am- 
bipolar diffusion is to gravitationally redistribute mass and magnetic flux within 
the interior magnetic flux tubes of a cloud, thereby initiating the formation and 
collapse of cores within magnetically supported envelopes. He further argued 
(Mouschovias 1979, 1982) that core formation will occur on the ambipolar diffu- 
sion timescale, r AD ~ Tg/r n j, where tq is the free-fall time and r n \ is the neutral- 
ion collision time in a cloud's inner flux tubes. Early studies (Mouschovias 
1979; Nakano 1979; Lizano & Shu 1989) focused on quasistatic evolution and 
the approach to core formation. 

The general theoretical scenario sketched above was verified by numeri- 
cal simulations of the evolution of two-dimensional axisymmetric (Fiedler & 
Mouschovias 1993) and flattened, thin-disk (Ciolek & Mouschovias 1994, 1995 
[CM94, CM95]; Basu & Mouschovias 1994, 1995a, b [BM94, BM95a, b]; Ciolek 
& Basu 2000 [CB00]), self-gravitating, isothermal, magnetic molecular cloud 
models. Following the evolution of the central flux tubes of a model cloud 
from the initial quasistatic contraction (due to ambipolar diffusion) to the later 
phase of dynamical collapse, from initial densities n niC o ~ 10 3 cm -3 to final 
densities n n<c ~ 10 11 cm -3 , these model simulations found that supercritical 
[(M/&b)cotc > (M/&B)crit\ cores did indeed form and separate within mas- 
sive, subcritical [(M/<E>£) cnv < (M/$>B)crit\ envelopes. The core masses found 
in these models were typically in the range ~ 3 — 30 M Q , with characteristic 
radii in the range 0.1 — 0.3 pc, depending on the choice of model parameters 
such as the initial central mass-to-flux ratio n c o (in units of the critical value for 
collapse) and the initial central degree of ionization xi iC o (= «i,co/ n n,cO> where 
n; iC o is the ion density). Mean core densities (n n ) were ~ 10 4 cm~ 3 . A magnetic 
disk-model was presented by Crutcher et al. (1994) to explain the evolution of 
the Barnard 1 (Bl) molecular cloud. Model predictions for the Bl core were 
found to be in excellent agreement with observations; predicted values of the 
core mass, and mean density and magnetic field strength were within 10% of 
the observed values. Further discussion of quantitative comparison of ambipolar 
diffusion models with observations can be found in § 3.1 below, as well as the 
paper by Bacmann in this volume. 

The disk- model clouds also included other interesting physical effects. Dy- 
namical effects of dust grains (charged and neutral) were studied in detail by 
Ciolek & Mouschovias (1993, 1994, 1995). Enhanced friction due to collisions 
with grains increases the coupling of neutral particles to the magnetic field, and 
can significantly lengthen the time required to form a core. Another interest- 
ing result from these models is that ambipolar diffusion reduces the dust-to-gas 
ratio in cores (also discussed in Ciolek & Mouschovias 1996), which, in turn, 
affects the ion chemistry in contracting cores, since ion and electron abundances 
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in the gas phase are dependent on the rate of capture on grains (an analytical 
derivation is provided by Ciolek & Mouschovias 1998). Rotation and magnetic 
braking were incorporated by Basu & Mouschovias (1994, 1995a, b). These 
studies showed that magnetic braking is generally very efficient (occurring on a 
timescale much shorter than r AD ) , allowing a core to form and collapse with very 
little hindrance by centrifugal forces; the angular momentum per unit mass J/M 
in these models was found to consistent with observations of rotating cores and 
estimates of likely values in the early solar nebula. The effect of the interstel- 
lar ultraviolet radiation field was investigated by Ciolek & Mouschovias (1995). 
They found that the increased ionization due to an external UV field greatly 
decreased the efficiency of ambipolar diffusion in cloud envelopes, significantly 
reducing the mass infall rate beyond the boundary of a supercritical core. 

Further details of the models can be found in the studies cited above, as 
well as the more recent reviews mentioned in § 1. Some more discussion of the 
L1544 model of CBOO is provided in § 3.1 below. 

3. Controversies, Unresolved Questions, and Work in Progress 
3.1. Extended Infall in Cores 

The protostellar core within the LI 544 cloud, located in Taurus, has received 
considerable attention in recent years by various workers. One such study was 
by Tafalla et al. (1998), who determined that the neutral gas in the core was 
infalling at a speed |t> n | ~ 0.1 km s _1 on a scale ~ 0.1 pc. A later study by 
Williams et al. (1999) used N2H + lines to find the infall speeds of the ionized 
gas component, and found |t>;| ~ 0.08 km s _1 on smaller scales ~ 0.02 pc. 
Both papers contrasted their findings with the ambipolar diffusion models of 
CM95 and BM94, and suggested that a discrepancy existed because certain 
ambipolar diffusion models did not exhibit such large infall speeds at extended 
radii. However, as noted by CBOO, many of the features observed in L1544 
(such as the ion infall speed and density structure on smaller scales) were in fact 
consistent with the previously published generic ambipolar diffusion models, if 
one considers the observed region in L1544 to lie within a supercritical core. Even 
in the absence of a specific model, magnetic restoring forces must be important 
in the evolution of L1544, since the magnitude of the inferred infall speeds 
are less than characteristic speeds such as the isothermal sound speed C in 
the gas (~ 0.2 km s _1 for this cloud) and the free-fall speed at these radii 
(S> C). And, as the linewidths in the L1544 core are very narrow (Williams 
et al. 1999), "turbulent" support cannot be invoked to explain the evidently 
retarded collapse. 

CBOO presented an ambipolar diffusion model specifically designed to study 
the evolution of the L1544 core. They modeled L1544 as a disk-like cloud inclined 
at an angle 6 = 16° to the plane of the sky; the resulting model evolution 
produced a core with radius 0.3 pc, which is indeed greater than the spatial scale 
studied by Tafalla et al. (1998) and Williams et al. (1999). Figure 1 (left) shows 
the spatial profile of the density at ten different times tj (see caption). Also 
shown (right) is the infall speed of the neutral gas v n (solid lines) and the ions V{ 
(dashed lines) at the same times tj. Qualitatively the overall behavior is similar 
to that seen in the earlier published ambipolar diffusion models, including the 
flattening of the density profiles at small radii, which is very similar to the density 
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profiles seen in sub-mm dust continuum studies of prestellar cores (e.g., Andre, 
Ward-Thompson, & Motte 1996). In Table 1, CBOO model predictions at times 
£2 and £3 are compared with various observed quantities (such as the mass M, 
density n n , radial neutral and ion speeds u n , r and Vi >r , and the mean line-of-sight 
magnetic field strength {B)\ os ) at several different radial locations. Given the 
likely observational uncertainties, we note that all of the model predictions are in 
reasonable agreement with the actual measured values — including the extended 
infall observations and the measurement of {B)\ os by Crutcher & Troland (2000), 
made subsequent to our prediction. [] This suggests that the current state t of 
L1544 is described by the theoretical model for £2 < £ < £3. 

The CBOO L1544 model is thus able to simultaneously explain the obser- 
vations of several disparate observations at many different radial locations. We 
also note that the magnetic models were fully developed prior to 1995, and were 
applied, without modification, to the recent observations of L1544 in a straight- 
forward fashion. That a single magnetic model can explain the observations 
with (e.g., L1544) and without extended infall through variations of input pa- 
rameters within the observationally allowable range, is particularly appealing, 
and in line with our discussion in § 1 of what attributes (simplicity, elegance, 
self-consistency, and predictive power) a proper physical theory should have. 

TABLE 1 
Physical Quantities in the L1544 Core 







CBOO: £ 2 


Observed 


CBOO: £3 


r ~ 0.14 pc : 


a 








M 




7.1 M 


8 M 


7.7 M 


|^n,r | 




0.09 km s" 1 
0.06 km s" 1 


0.1 km s _1 


0.11 km s" 1 
0.07 km s" 1 


\ v i,r\ 




r ~ 0.06 pc : 


b 








(-B)los 




13.4 fiG 


11 fiG 


15.5 fiG 


r ~ 0.02 pc : 


c 


0.07 km s" 1 




0.14 km s _1 


|^n,r | 




\ v i,r\ 




0.05 km s" 1 


0.08 km s" 1 


0.11 km s" 1 


n n 




2.3 x 10 5 cnr 3 


4 x 10 5 cm -3 


3.3 x 10 5 cnr 3 


r ~ 1.2 x 10 3 AU : ' 


l 






n n 




4.1 x 10 5 cm~ 3 


7 x 10 5 cm" 3 


2.2 x 10 6 cm" 3 



a Observed values taken from Tafalla et al. (1998). Observed values taken 
from Williams et al. (1999). c Observed value taken from Crutcher & Troland 
(2000), measured after the prediction made by CBOO. d Observed value taken 
from Ward-Thompson et al. (1999) and Bacmann et al. (2000). 



J The values for the drift speed v D (— Vi — v n ) which lie in the range 0.02 km s 1 to 0.04 km s 1 
for the L1544 model, are also in agreement with the limits deduced by Benson et al. (1998), 
who inferred v D <. 0.03 km s _1 from velocity differences of neutral molecules (CCS and C3H2) 
and NH3 + in sixty dense cores. 
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Figure 1. Spatial profiles of physical quantities in the LI 544 model 
cloud of CB00, at ten different times tj (j=0, 1, 2,.. .,10), at which 
the central density at time tj is a factor lCP greater than the initial 
central density. These times are, respectively, 0, 2.27, 2.60, 2.66, 2.680, 
2.684, 2.685, 2.6856, 2.68574, 2.68577, and 2.68578 Myr. An asterisk 
on a curve (present for tj > to)) locates the instantaneous radius of the 
critical magnetic flux tube. Left: density. Right: infall speed of the 
neutrals (solid curves) and the ions (dashed curves). 
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Figure 2. Initial flux constant C(/i c o) = t COTe /T 9c0 for several model 
clouds with k = 0.5, 0.6, 0.75, and 1.0, respectively. Also shown are the 
values found from previously published numerical simulations: model A 
of CM94 (located by the center of ®), model 2 of BM94 (©), model 5 of 
BM95b (•), and the L1544 model cloud of CB00 (o). Above fi c0 £ 0.5, 
C(a*co) * s significantly less than unity. 
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3.2. Timescale for Star Formation 

From the observed infall data discussed above for L1544, one may deduce a 
kinematical timescale ~ r/|u n | ~ 1 Myr. Recent statistical studies of protostellar 
cores and young stellar objects (Lee & Myers 1999; Jijina, Myers, & Adams 
1999) have also estimated a similar evolutionary timescale for star formation. 
These studies also suggested that the observed timescales were in disagreement 
with earlier published ambipolar diffusion models (e.g., CM95; BM94), which 
had cores forming on a timescale t cor e — r $ c0 ~ lOrff, where r $ c0 = r AD c0 /2 = 
(r|/2r n i) c o is the initial central flux-loss timescale, as proposed early on by 
Mouschovias (see § 1). 

However, the relation t corc ~ r $ c0 is strictly valid only for highly subcritical 
model clouds, i.e., clouds with ^ c q < 1. In fact, ambipolar diffusion models 
with less subcritical initial mass-to-flux ratios have t cor e *C t # . For instance, 
the L1544 model of CBOO had // c o = 0.8 and t core = 0.2t $c0 ; numerically this 
was equal to 1.3 Myr, similar to the values from the statistical studies cited 
above. Hence, t core is dependent on /x c q. Ciolek & Basu (2001; [CB01]) revisited 
the problem of core formation by ambipolar diffusion, and derived an analytical 
expression for t core valid for all subcritical clouds (fi c o < 1). Analyzing the evo- 
lution during the epoch of core formation, they found that t core = C(// c o)t #c0 , 
where C(/i c o) is the initial flux constant of a cloud, and is a measure of "how 
far away" a cloud is from the critical mass-to- flux ratio (i.e., fi c = 1 in dimen- 
sionless units). For clouds in which neutral-ion collisions are predominant, and 
approximating the relation between ions and neutrals in the central flux tubes 
of a cloud by n i)C ex n^ c (0.5 £ k £ 1, Ciolek & Mouschovias 1998; [CM98]), 



C (M ,r°^ forfc<l, (1) 

for Jfe = 1 (2) 



(a derivation is given in § 2.1 of CB01). 

Figure 2 shows C(// c o) for several model clouds with different values of k 
(see caption). Also shown are the values of C(// c o) = *core A" $c0 f° r previously 
published ambipolar diffusion models; the results from the numerical simula- 
tions are in excellent agreement with the analytical expressions (|l|) and (g). For 
most realistic values of k, clouds with ^ c o <C 1 have C(/x c o) ~ 1; an d, therefore, 
*core ~ 7" Sc0 > as Mouschovias (1979) originally argued. However, for /i c o ^ 0.5, 
C(Mco) "C 1, with C{fjL c o) ~ (1 — ^ c o) for 1 — /i c o <S 1. This region of param- 
eter space is especially of interest, since, as noted in § 1, Crutcher (1999) has 
measured mass-to-flux ratios in several clouds consistent with /i c ^ 0.5. While 
it is likely that Crutcher's data set contains clouds that are evolving towards or 
have recently exceeded the critical state, and are therefore beginning to collapse 
(which is what we suggest has occurred in L1544; see § 3.1 above), it may also 
contain typical values of fj, C Q. This means that t core will be <C r # c0 for clouds 



1 - u {1 ~ k) 

1 McO 


2(1 -it) 

V^cO/ 
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with these particular values o//x c o. Numerically, r #c0 ~ 7 (xi. c o/10~ 7 ) Myr, f] 
where x- 1)C q = ni jC o/n njC o has been normalized to values typical of the interior 
of molecular clouds for 10 3 cm~ 3 <J n n>c <^ 10 4 cnr 3 (e.g, CM98). It follows 
then that t COTe = C(^ c o)r #c0 <; 4 Myr for /i c o > 0.5. Thus, the recent statistical 
analyses suggesting that star formation occurs in some regions on timescales ~ 
1 Myr may be explained by observations of subcritical clouds and ambipolar dif- 
fusion models with central mass-to-flux ratios slightly below the critical value for 
collapse. Further discussion of the consistency and implications of ambipolar 
diffusion models with estimated timescales for protostellar evolution, including 
the effect of dust grains, is presented in CB01. 

3.3. Core "Edges" 

Bacmann et al. (2000) have recently conducted an extensive ISOCAM absorp- 
tion survey of prestellar cores, including L1544, and made quantitative compari- 
son of their results to theoretical models of core formation and collapse. For the 
most part, they found that the ambipolar diffusion models of CM95 and BM94 
provided the best model fits to their data in the interiors of cores. However, in 
some cores they found that on larger scales the cores effectively had "edges", 
where the density decreased rapidly with increasing radius — much more so than 
predicted by the theoretical models. For L1544, r e d ge ~ 0.14 pc. fj Bacmann et 
al. suggested that the presence of these edges indicates that the current theoret- 
ical models neglect some important physics on larger scales, such as turbulentR 
support, which has not yet been incorporated into the ambipolar diffusion sim- 
ulations. However, if the turbulence were due to hydromagnetic waves (HMW), 
as advocated by Mouschovias (1987) and Mouschovias & Psaltis (1995), then 
it is possible that the reason the non-turbulent CB00 model matches observa- 
tions for r ^ r e d gc in LI 544 is due to the decay of the turbulence for scales less 
than the Alfven minimum lengthscale Aa = ttuaTiu (Kulsrud & Pearce 1969); 
va = -B/(47T/9 n ) 1 / 2 is the Alfven speed in the neutral gas. As first proposed by 
Mouschovias (1987), ambipolar diffusion will damp HMW for scales < Aa, result- 
ing in a loss of support interior to this region when the HMW turbulence decays 
away. Formation and collapse of cores can then proceed via the 'standard' self- 
initiated ambipolar diffusion models, described earlier. Numerically, the expres- 
sion for the damping length Aa = 0.41(_B/10/xG)(10 3 cm~ 3 /n n ) 3 / 2 (10~ 7 /xi) pc, 
which, when evaluated for the initial parameters of the CB00 L1544 model, is 
~ 0.1 pc. Hence, if HMW were to be important in L1544, we would expect 
that a non-turbulent ambipolar diffusion model would deviate from the actual 
evolution on scales ^ Aa ~ ^edgej as observed. Moreover, if decay of HMW has 
occurred for scales < Aa, the linewidths should be narrow in the core, which is 
again observed in LI 544. 



2 Note that x- uc o can be as small as 10 8 (CM98) so that a low ii, c o is another means of reducing 
the flux-loss timescale (e.g., see models 8 and 9 of BM95a), and, hence, t CO rc 

This value is after deprojection. The measured value for r d gc quoted by Bacmann et al. (2000) 
is ~ 0.04 pc. 

4 In this paper, we use the term "turbulent" quite generally, to refer to all forms of non-thermal 
motions, whether an ensemble of wave motions or a true turbulent cascade. 
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3.4. Protostellar Accretion Rates 

The later stages of core collapse has also been investigated. Basu (1997) stud- 
ied semi-analytically the approach of collapsing cores to higher densities than 
studied in the earlier numerical simulations. He suggested that ambipolar diffu- 
sion would continue to affect the increasingly dynamical infall even during the 
higher density stages of evolution, and also that the accretion rate M would 
become significantly larger than the well-known, time-independent expression 
M ~ C 3 /G = 1.6 (T/10K) 3 / 2 M Myr" 1 found in self-similar collapse models 
of a singular isothermal sphere (Shu 1977). This was confirmed by Ciolek & 
Konigl (1998; [CK98]), who followed the collapse of a magnetic protostellar core 
(including the effect of ambipolar diffusion) up to and beyond the formation of a 
finite mass object at r = 0, which they dubbed "point-mass formation" (PMF). 
They followed the evolution in their models up to the time at which a 1 M & ob- 
ject had accumulated at the center of a core. The accretion rate in their typical 
model reached a maximum value of 9.4 M Q Myr~ near the onset of PMF, and 
then declined to approximately half that value by the time a 1 Mq protostar 
had formed. During the later stages of the evolution of their model, ambipolar 
diffusion affected the accretion rate through the formation of a hydromagnetic 
shock at radii r ^ 100 AU; the relevant physics of the inflow and accretion of the 
numerical model was also captured in a companion self-similar collapse study 
presented by Contopoulos, Ciolek, & Konigl (1998). 

The duration of the CK98 simulations spanned the evolution of a proto- 
star from a class to a class I object. They suggested that the time-dependent 
accretion occurring in their models was in accordance with the analysis of Bon- 
temps et al. (1996), who proposed that a decline in the observed momentum 
flux of CO outflows from young stellar objects with decreasing envelope mass 
indicated a time-dependent M during this transitional phase of a protostar's 
life. CK98 also compared their results to the measurements of the gas in the 
infalling disks of HL Tau (Hayashi, Ohashi, & Miyama 1993) and L1551-IRS5 
(Ohashi et al. 1996). The late-time accretion rate and gas mass of the CK98 
model at r ;> 500 AU were consistent with observed values for both objects. 

3.5. Magnetic Field Strengths and Alignment 

A remarkable feature of the magnetic field strength measurements compiled 
by Crutcher (1999) is that they are consistent with the intrinsic mass-to-flux 
ratios (in units of the critical value) of all cores being very close to unity (see 

1/2 

earlier discussion in § 1). Additionally, the observed correlation B oc p n , 
which implies some degree of flattening along the mean magnetic field direction, 
can even be improved if one accounts for internal support due to non-thermal 
motions (Basu 2000). The above two features of the measurements implies that 
cores have dynamically important magnetic fields and are characterized by force 
balance along magnetic field lines. 

Recent submillimeter polarization maps of dense regions of molecular clouds 
(Matthews & Wilson 2000; Ward-Thompson et al. 2000), in which the polariza- 
tion vectors truly sample the dense gas, reveal that the magnetic field direction 
is invariably either parallel to the apparent minor axis of the objects or has 
a small offset angle \I/ relative to it. Ward-Thompson et al. (2000) measured 
\l/ ~ 30° for three cores (including LI 544) and wondered how strong magnetic 
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Figure 3. Simulated contours (solid lines) and mean polarization 

direction (dashed lines) for a triaxial body with axial ratios £ = c/a = 
0.3 and £ = b/a = 0.6 seen from three sets of viewing angles (9, (ft). 
The The x' and y' axes are in the plane of the sky and are chosen to 
lie along the projected major and minor axes, and a' is the apparent 
semi-major axis. The dashed lines lie along the projected direction 
of the intrinsic z-axis. a) (#,</>) = (tt/2,0), yielding an apparent axis 
ratio q = 0.3 and offset angle ^ = 0°. b) (9,4>) = (vr/4, 7r/4), yielding 
q = 0.57 and * = 28°. c) (9, <f>) = (vr/6,0), yielding q = 0.68 and 
^ = 90°. 



fields were compatible with any nonzero *$>. The resolution to this issue comes 
from the fact that the axisymmetric assumption in the magnetic models must be 
relaxed, even while the assumption of a dynamically important magnetic field 
is maintained. As pointed out by Basu (2000), any non-axisymmetry in the 
core shape (which makes them triaxial bodies) results in a finite probability of 
viewing a nonzero ^f. Figure 3 illustrates how a single object can appear to 
have any possible value of \&, depending on the viewing angles (9, (/)) defined on 
a spherical viewing surface around the triaxial body. However, if the shortest 
axis of the body is along the mean magnetic field direction, there is an expected 
correlation towards ^ = for a large sample of measurements (Basu 2000). 
Hence, large data sets can easily distinguish between this scenario and that ad- 
vocated by Ballesteros-Paredes, Vazquez-Semadeni, &: Scalo (1999), who use a 
highly turbulent model (with a weak magnetic field) to claim there should be 
no correlation between cloud elongation and magnetic field direction. 

The above explanation for the offset *$> = 29° observed in L1544 by Ward- 
Thompson et al. (2000) means that the axisymmetric assumption employed by 
CB00 needs to be relaxed, if only slightly. The correspondence of the Crutcher 
& Troland (2000) field strength measurement with the prediction of the CB00 
model, and the overwhelming evidence for flattening along the mean field direc- 

1/2 

tion in a variety of cores (B oc p Q ' ) implies that the oblate geometry of the 
CB00 model is a reasonable approximation. 

3.6. Turbulence and Cluster Formation 

Individual regions which form stars likely must damp their internal turbulent 
motions and also overcome the support due to the mean magnetic field. Since 
these motions may damp preferentially at high density or column density where 
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the coupling of charged particles is reduced, and the mean magnetic field ef- 
fectively introduces a threshold column density for collapse (§ 2), a fruitful 
approach to understanding cluster formation may be to study the properties of 
peaks above a critical threshold in a turbulent column density field. Basu & 
Pudritz (2000) have studied the properties of a Gaussian field of column density 
fluctuations and how it relates to the properties of observed clusters. As they 
point out, star formation in terms of thresholds, rather than say, a preferred 
scale of fragmentation, is qualitatively in better agreement with observations. 
The threshold idea leads to the result that regions of high mean column den- 
sity (close to a critical threshold) are likely to form stars (in regions exceeding 
the threshold) that are more closely spaced and typically more massive than in 
regions of low mean column density (far from the threshold) . In linear fragmen- 
tation theory, regions of high column density yield a smaller fragmentation scale 
as well, but also favor stars of relatively lower masses. 

Finally, the pressure due to an external column of (turbulent) matter may 
play a role in limiting the size of infall regions and creating core "edges" at 
smaller radii than in models without external turbulent motions (see related 
discussion in § 3.3). Despite objections to the concept of a confining pressure 
raised by Ballesteros-Paredes et al. (1999) on the basis of a single simulation, we 
note that an external pressure must exist simply due to the weight of an external 
cloud or cloud complex. While turbulent motions obviously exert kinetic effects 
on some scale, neglecting external pressure is as absurd as neglecting the Earth's 
atmospheric pressure on a windy day! 

4. Summary 

In this paper we have described several aspects of the theory of star formation 
in magnetic interstellar clouds, first formulated in the 1970's and developed 
further in the past decade. This single theoretical framework is consistent with 
(and can predict) various seemingly disparate facts such as the inefficiency of 
star formation, the structure of regions of core formation and collapse, magnetic 
field strengths and mass-to-flux ratios, transport of angular momentum, time- 
dependent mass accretion rates, magnetic field orientations in clouds, etc. In 
fact, specific models for core formation and evolution in two molecular clouds 
(Barnard 1 and L1544) have already been presented, and found to be in good 
agreement with observations for these two clouds. We have also discussed how 
recent data related to extended infall, timescales for star formation, and the 
alignment of magnetic fields in cores, are compatible with the results of magnetic, 
ambipolar diffusion models. Our approach is to make quantitative predictions 
and, where possible, clarify differences with other competing scenarios. For 
instance, sub-mm polarimetric observations indicating some degree of alignment 
of magnetic fields and core axes strongly favor magnetic, ambipolar diffusion 
models, rather than highly turbulent star formation models advocated by some 
other workers, which predict no such correlation. It is through such rigorous 
testing that a scientifically valid theory of star formation will ultimately be 
determined. 
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